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Abstract 

We present data-driven methods for the full reconstruction of jets in heavy ion collisions, for inclusive and coincidence jet 
measurements at both RHIC and LHC. The complex structure of heavy ion events generates a large background of combinatorial 
jets, and smears the measured energy of the true hard jet signal. Techniques to correct for these background effects can induce 
biases in the reported jet distributions, which must be well controlled for accurate measurement of jet quenching. Using model 
studies, we evaluate the proposed methods for measuring jet distributions accurately while minimizing the fragmentation bias of 
the measured population. 
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1. Introduction 

The interaction of QCD jets with the Quark-Gluon Plasma (QGP) generated in high energy nuclear collisions ("jet 
quenching") provides unique and penetrating probes of the QGP. Jet quenching has to date been measured primarily 
via the suppression of high pj hadrons and their correlations [1]. Full jet reconstruction (2l[3) may enable the study 
of jet quenching at the partonic level, without the added complexity of hadronization. However, jet reconstruction 
in heavy ion events is a challenging task, due to the large population of combinatorial background jets generated 
by the jet reconstruction algorithm from random recombination of particles not correlated via single hard scattering, 
and to the distortion of the measured energy of true jets by background fluctuations. Experimental control of these 
effects may be achieved by pruning the soft component of events prior to jet finding, and requiring jet candidates to 
contain a hard fragmentation component. However, for quenched jets, these analysis techniques may also bias the jet 
population and the reported jet energy. For accurate quenching measurements it is necessary to assess such biases 
and to minimize their effects. In these proceedings we present data-driven methods for minimally -biased inclusive 
and coincidence jet measurements in heavy ion collisions, at both RHIC and LHC. We base our approach on previous 
developments |H[5][6]|. Background fluctuations are "unfolded" using an iterative technique incorporating Bayes's 
Theorem [7]. We employ a model event generator to evaluate these methods, and give prescriptions to minimize 
fragmentation bias in real data analysis. This approach is applicable to jet measurements by the STAR experiment at 
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RHIC and ALICE at the LHC, which are able to reconstruct a large fraction of all jet fragments via precise charged 
particle tracking and electro-magnetic calorimetry. 

2. Toy Model and Jet Reconstruction 

Our aim is to clarify the generic effects underlying the largest systematic uncertainties in heavy ion jet measure- 
ments. For such a parametric study, it is not necessary to model heavy ion events and experimental response in 
detail; indeed, the complexity of such an approach may obscure the important generic effects. We therefore employ 
a simplified Toy Model event generator, that nevertheless captures the main features of heavy ion jet reconstruction 
seen in data (UEIE). The Toy Model randomly generates massless particles with uniform distribution within \rj\ < 1 
and < (fi < 2n. The transverse momentum distribution of the soft background is a Boltzmann distribution, with 
(p T ) = 500 MeV at RHIC and (p T ) = 700 MeV at the LHC. The choice of jet distribution depends on the observable: 
for inclusive jets we utilize the inclusive spectrum measured in p+p collisions at the same scaled by Taa, while 
for the coincidence measurement we utilize the semi-inclusive distribution calculated by PYTHIA for p+p collisions, 
without additional scaling. Jets are fragmented using the PYTHIA fragmentation routines, or are not fragmented at 
all; i.e. a single particle carries the total jet momentum ("SP" fragmentation). The total multiplicity per event corre- 
sponds to the average multiplicity of all particles (charged plus neutral) in 0-5% central Au+Au collisions at RHIC or 
0-5% central Pb+Pb collisions at the LHC. 

Millions of Toy Model event are generated, and are analyzed with the same algorithms used for real data. Jets are 
reconstructed using the Anti-^^ algorithm [ 8 ] with resolution parameter R - 0.4 and energy recombination scheme, 
and with median background density estimate p as described in (9) . We utilize the FASTJET implementation with 
default settings ifTOl . Jet reconstruction incorporates all particles with p T > 0.2 GeV/c and \r]\ < 1.0, with full 
azimuthal coverage. Accepted jets have \q\ < 0.6 and area Aj et > 0.4 (4). The measured jet energy p™ c is corrected on 
an event- wise basis via (9) : 

p (™> = p ™- p . Ajet , a) 

where p is the event-wise estimate of the background density (energy per unit area) and (corr) indicates event-wise 
correction for background. Local variations of background density relative to p will generate large distortions in 
the distribution of p^ orr \ which must be corrected via unfolding. We use iterative unfolding incorporating Bayes's 
Theorem (7), with the response matrix corresponding to Spj, the distribution of fluctuations in jet response measured 
using data [fflEl. In this case, the "data" are Toy Model events. 

The results presented here are for SP fragmentation, in order to isolate the biases in the jet measurement due to 
suppression of the combinatorial background component and unfolding of background fluctuations. The validity of 
the SP approach is based on the observed insensitivity of the response of the Anti-&r algorithm in heavy ion events 
to the detailed pattern of fragmentation of particles measured within the jet cone |H|5)- Figure [T] shows the generated 
particle distribution for SP fragmentation, together with the reconstructed (/?^ orr) ) jet spectrum, for central Au+Au 
collisions at RHIC. 
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Figure 1. Toy Model distributions for central Au+Au collisions at RHIC. Left: particle distribution with SP fragmentation. Right: particle ("Truth") 
and reconstructed jet spectra. 



10 3 





RHIC Kinematics 

Boltzmann Distribution 






Jet Disitribution 
Particle Disitribution 




































I 4 6 £ 


? '10 



10 3 
10 2 
10 

1 

10" 1 
10 2 
10"" 
10" 
10" 
10- 
10" 
10" 



_ • Truth 

Measured 



P T (GeV/c) 



2 



G. O. V. de Barros, Bo Fenton-Olsen, Peter Jacobs, Mateusz Ploskon / Nuclear Physics A 00 (2012) 7-^| 



3 



10 3 



RHIC Kinematics 


— Truth 


- u-o~/o central collisions 

N evBn , s = 10M 

Anti-k T R=0.4 
" A Jet > 0.4 sr 




Measured 




. Unfolded 
Unfolded 


N. =1 
N' ter = 2 

iter 




0.2 GeV/c " T'°°<, 























































20 30 40 
(GeV/c) 



LHC Kinematics 



10 2 -o- 5% central collisions 

10Av e n,s=1M 

1 Anti-k T R=0.4 
n -1 ^jet * 4 sr 

10 =| , eading >02 

10 ? 
10" 3 _ 
10" 4 = 
10" 5 - 

10" 6 _ 

10" 7 = o° 

10 -100 -50 



Truth 
Measured 
Unfolded, N. 
Unfolded, N* e 



50 100 

p™ rr (GeV/c) 



10 3 
10 2 
10 
1 

10" 1 
10" 2 
10" 3 , 

; 10 3 
j 10 2 
j 10 

! 1 

10" 1 

10~ 2 
10" 3 f 



zRHIC Kinematics 
0-5% central collisions 

^ N e v B „ts = 10M 

Anti-k T R=0.4 
, A, e , > 0.4 sr 

Reading > g 2 GeV / C 



. Unfolded, N jter = 1 
. Unfolded, N. ter = 2 



30 35 40 
P T (GeV/c) 



|_HC Kinematics | 
-0-5% central collisions 

^events = 1M 

:Anti-k T R=0.4 
A, et > 0.4 sr 

Reading > Q 2 GeV/( . 



. Unfolded, N. (c 
. Unfolded, N. 



100 120 
P T (GeV/c) 



Figure 2. Measured, truth and unfolded distributions for fully inclusive jets. Upper: RHIC; lower: LHC. Right plots show ratio Unfolded/Truth. 



3. Inclusive Jet Spectrum Measurement 

We first consider measurement of the fully inclusive jet spectrum. Figure [2] shows the Truth, measured, and 
unfolded jet distributions for RHIC (upper panels) and LHC (bottom panels). The distributions resulting from the 
first two Bayesean unfolding iterations disagree strongly with the Truth distribution, and successive iterations do not 
reduce the discrepancy. This effect is due to the presence of an overwhelming population of combinatorial background 
jets which do not have an underlying physical distribution, but which the unfolding algorithm cannot distinguish from 
true jets. Figure [3] shows the same calculation, but with the distribution to be unfolded restricted to jets containing 
a leading particle /^ ading > 4.0 GeV/c at RHIC and /^ ading > 10.0 GeV/c at LHC. This cut strongly suppresses the 
combinatorial jet background distribution, and the resulting Unfolded distribution corresponds to the Truth distribution 
within 10% for all pj. 
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Figure 3. Measured, truth and unfolded distributions for inclusive, jets with leading particle above a given threshold. RHIC (upper plots) and LHC 
(bottom plots). 
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The calculations presented here employ simplified SP fragmentation, though calculations employing PYTHIA 
fragmentation show similar behavior as a function of the /?^ admg threshold. The behavior in data will depend on the 
detailed interplay of (quenched) jet fragmentation and the background hadron distribution, including effects of flow, 
which are not included in the Toy Model. 

It is apparent that accurate unfolding of background fluctuations from the inclusive spectrum requires suppression 
of the combinatorial background, via imposition of a fragmentation bias (see also ifTTIl ). However, this approach also 
biases the resulting true jet population, which may be significant for quenching measurements. The lowest threshold 
value for p l ™ dmg that enables stable unfolding of the background fluctuations is clearly preferred, but quantitative 
assessment of the remaining biases in such an approach are beyond the scope of this model study. 

Given effective suppression of the combinatorial jet population, the precision of the inclusive measurement is 
then limited by the precision with which the unfolding response matrix is known. The response matrix includes 
contributions both from detector effects and from the fluctuating event background, with the latter generally dominant 
in heavy ion measurements. Data-driven techniques ElOEl are able to measure the fluctuation distribution Spj over 
several decades in magnitude. 



4. Hadron+ Jet Coincidence Measurement 

We next turn to a coincidence measurement, the semi-inclusive rate of jets recoiling against a high hadron 
trigger ("h+jet"). The acceptance for the recoil jet is 

l^hadron ~ <£jet ~ 7l\ < 7T/4. (2) 

The basic observable is the pj distribution of the number of observed recoil jets, normalized by the number of triggers. 
We prefer a hadron to a jet trigger, since hadrons can be accurately measured in heavy ion collisions without having 
to account for complex background effects, and their (suppressed) production is well-understood in jet quenching 
models. In addition, model studies suggest that a high p T hadron trigger imposes a significant "surface" bias on the 
measured population, corresponding to maximum the path length in matter for recoiling jets, whereas a jet trigger 
does not (T2H3. 
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Figure 4. h+jet measurements at RHIC (yfs = 200 GeV) : semi-inclusive jet distributions. Error bars are from event generation and correspond 
to very large integrated luminosity. Upper left: PYTHIA calculation for p+p collisions with various intervals of trigger hadron pj. Remaining 
panels are Toy Model calculations for central Au+Au collisions. Upper right: p^ on ^ distribution for two intervals of trigger hadron pj; lower left: 
difference distribution and results of unfolding; lower right: ratio Unfolded/Truth. 

Fig. |4j upper left panel, shows a PYTHIA calculation of the semi-inclusive h+jet distribution for p+p collisions at 
RHIC = 200 GeV, for three different intervals of hadron trigger /? T . Note that the recoil jet distribution depends 
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strongly on hadron trigger /? T , and the probability to observe a recoil jet in the acceptance with p T above the hadron 
trigger threshold is at most a few percent. Fig. [4j upper right, shows the equivalent recoil jet distribution for central 
Au+Au collisions, for two intervals of trigger hadron pj. Unfolding of background fluctuations for these individual 
distributions fails in a way similar to unfolding of fully inclusive distributions (Fig. [2]), and for the same reason: the 
recoil region is dominated by combinatorial background jets, even for a high pj trigger, due to the low probability for 
true coindiences. In order to overcome this background effect without imposing a fragmentation bias on the recoil jet 
population, we note that the combinatorial jet distribution is, by definition, independent of trigger /? T . This invariance 
is seen in the upper right panel, which compares the recoil p^ ou) distribution for two different trigger pj ranges. The 
region p^ on) < 0, where the combinatorial background component may be expected to dominate, is very similar for 
the two trigger intervals, whereas the region for large positive p^ on) > exhibits strong correlation with trigger p T . 
This invariance may not be precise in real data analysis, due to biases in reaction plane orientation and centrality 
imposed by the requirement of a high p T trigger hadron. Such biases are not present in Toy Model events. Their 
effects in real data analysis can however be minimized, as described below. 

We exploit this invariance by defining a new observable: the difference between the two distributions in Fig. [4] 
upper right panel, which represents the evolution of the recoil distribution with trigger /? T . The hadron trigger /? T for 
both ranges should be chosen to be high enough that the probability per event for such a trigger is low, so that a hard 
recoil jet most likely originates from the same hard interaction as the trigger. In addition, for sufficiently high pj 
hadron triggers, the reaction plane and centrality biases of the trigger have weak, if any, dependence on trigger p T , 
and such biases will affect the combinatorial background similarly for the two trigger pj ranges. After accounting 
for the strict conservation of jet density Q, the combinatorial background jet component can be suppressed in the 
difference distribution in a purely data-driven way, at the per-mil level of precision or better. Fig. [4] lower left, shows 
the difference distribution ("Measured", red points), which represents the hard recoil jet distribution, but with energy 
still smeared by background fluctuations. The lower left panel also shows the result of unfolding of these fluctuations 
compared to the Truth distribution from PYTHIA p+p events, with the ratio Unfolded/Truth shown in lower right 
panel. The Truth distribution is recovered within ~ 10% precision over a broad kinematic range above the trigger 
threshold, with minimal cut on the jet constituents (pj > 200 MeV) and without the imposition of any fragmentation 
bias on the measured jet population. Similar results are obtained for h+jet measurements at the LHC, with appropriate 
choice of trigger p^. 
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